Proton Electro-Magnetic Form-Factors in the Instanton Liquid Model 
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We report the first non-perturbative calculation of proton electro-magnetic form-factors in the 
Random Instanton Liquid Model (RILM) and in the Interacting Instanton Liquid Model (IILM). 
By calculating the ratio of appropriate three-point to two-point functions, we divide out the cou- 
pling constants and compare our results directly to some integral of the form-factors. Using various 
parametrizations of the electric form-factor Ge{Q^) at large > 3.5 GeV^, where it is not yet 
measured, we compare those with expected theoretical dependence. We find from this comparison 
that some distributions of charge are clearly excluded (e.g. the same as of the magnetic moment, 
^Ge/Gm = 1, as well as the opposite scenario in which ^.Ge/Gm rapidly approach zero), restrict- 
ing possible behavior of the form-factor to rather narrow band. Furthermore, we found that our 
calculation of the nucleon form-factors is dominated by a configuration in which two out of three 
quarks interact with a singie-instanton, in spite of the fact that the evaluated three-point function 
has rather large distances ( ~ 1.2/m ) between points. We also estimate the size of the scalar 
di-quark and found it to be very small, comparable to the typical instanton size. 

PAC number(s): ll.lStk 



I. INTRODUCTION 

The non-perturbative sector of QCD is not yet theo- 
retically understood, and any direct comparison between 
high accuracy experimental data and theoretical predic- 
tions is very valuable. The electro-magnetic form-factors 
of hadrons are fundamental observables, depending on 
one scale Q^, which are related to the probability of 
an hadron to absorb a photon and yet remain the same 
hadron. Apart of telling us much about hadronic sizes 
and structure, and about the dynamics responsible for 
the ricombination of partons in the final state, form- 
factors are also supposed to teach us where the transition 
from the non-perturbative to perturbative regime takes 
place. 

Renewed interest to these issues is triggered by the 
new precision measurements pursued at the Jefferson 
Laboratory. Specifically, their new measurements (us- 
ing a polarization-transfer method) have yielded more 
accurate proton form- factors, showing for the first 
time that distributions of charge and magnetic mo- 
ment are not identical, ^Ge/Gm H- The ratio 
hGe{Q'^)/Gm{Q^) was measured up to Q^= 3A7 GeV^ 
and new data up to ~ 5.6 GeV^ should be available 
soon. These result are plotted at fig. |^ below. 

Theoretical arguments about this issues are rather 
controversial. Apparent agreement between existing 
data and pQCD power counting |], ig {F^{Q'^) ~ 
1/Q^,Gn ~ ^/Q*) have lead some to think that pertur- 
bative regime is reached (for a review of pQCD analysis 
of hadronic form- factors see, e.g. [Q). However, at least 
in the case of the pion form-factor we do know that the 
observed value of F.,r{Q'^) is different from its asymp- 
totics at infinite , roughly by factor two, and thus some 
adjustment of its value at higher than measured so 
far should take place. 

There are also numerous arguments which show that 
the origin of the form-factors at Q'^ < 2 GeV^ is non- 



perturbative. Historically, the first calculations of this 
kind for the proton form-factors were based on the Op- 
erator Product Expansion and QCD sum-rule analy- 
sis, supplemented by the assumption of vacuum dom- 
inance 1^,^,^. In this approach all information about 
the non-perturbative properties of the QCD vacuum are 
parametrized via few quark and gluon condensates. How- 
ever, this approach is only valid at small enough dis- 
tances, and in order to extract physical observables one 
is then forced to consider complete experimental spectral 
densities which includes contributions from multiple ex- 
cited states, not just pion or nucleon alone. To overcome 
this difficulty, a simple "pole-plus-continuum" (PPC) 
model for the spectral representation is often adopted. 
Unfortunately, that introduces significant model depen- 
dence in the analysis and reduces its predictive power. In 
practice, this approach is not always successful: for ex- 
ample the sum-rule calculation of the proton's two-point 
function involving a PPC ansatz agrees with the exper- 
imental data, only if a particular loffe current is used 
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FIG. 1. JLAB results for ^Ge/Gm as compared to 
the fits (E3) with 
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In general, hadronic data usually tell us about long- 
distance behavior, while field theoretic calculations have 
better predictive power in the opposite limit of small dis- 
tances. One issue any application of such method has to 
face is the existence and size of the "working window" 
of distances in which a quantitative comparison between 
the two can be made. To be successful, a theoretical ap- 
proach has to be able to predict the correlation functions 
at relatively large (~ 1/m) distances: then the unknown 
contribution from all excited states is automatically sup- 
pressed and only the pole corresponding to the lowest 
lying particle survives. Two illustrative examples of rela- 
tive strength of both components, for two-point isovector 
pseudoscalar current correlator ( "pion" ) and the nucleon 
(with current to be specified below) are shown in Fig.^ 
as a function of (Euclidean) distance I between the cur- 
rents. Here and below, all correlators are normalized to 
their free quark value (so unit value at small I means 
just free quark motion, a simple consequence of asymp- 
totic freedom). Frome those examples one can see, that 
the pole dominance is rapidly reached at some (channel- 
dependent) value of I, about 0.8 /m for the nucleon. 
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FIG. 2. Pole and continuum contributions for the pion 
isl (dashed, so = l.GGeV^) and nucleon (dot-dashed, 
2.5 GeV^) two-point correlation functions. 
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The most direct way of calculating ab initio such cor- 
relation functions is numerical lattice-based approach. 
However, quite successful results were also obtained from 
much simpler instanton-based models, see review |^ and 
also for recent example of quantitative predictions of 
the model, compared to the t lepton decay data. Closest 
to the present work is the successful calculation of the 
pion form-factor p^ , p^ using the same instanton-based 
model as used below. One important aspect of is the 
direct relation found there between the pion and instan- 
ton sizes. 

In the present work we try to answer the following 
questions. How to organize the calculation, what are the 
best combination of three and two-point correlators to 
be used? At what distances those should be calculated, 
in order to exclude contributions of the nucleon excited 
states? How do the results obtained correlate with exper- 



imental form-factors? Are they sensitive to the high-Q^ 
part which is not yet measured, and if so, can we make 
some predictions? Can the proton form-factors be de- 
scribed in terms of configurations of three valence quarks 
bound by instanton forces? Are these dominated by one 
single instanton, at least at some scale? If so, do all 
three quarks interact strongly with it, or only two, as ar- 
guments based on Pauli principle for zero modes would 
suggest? And finally, how does this approach correlate 
with phenomenological models of the nucleon, such as 
quark-diquark model? 

The paper is organized as follows. In section ||, we 
will present our results for the nucleon two-point and 
electro-magnetic three-point function, obtained in the 
Random Instanton Liquid Model (RILM), in the Inter- 
acting Instanton Liquid Model (IILM) and in the Single 
Instanton Approximation (SIA). Is section [H we show 
how such correlation functions can be evaluated from the 
phenomenological knowledge of proton mass and form- 
factors and we compare them with our theoretical predic- 
tions. In section we discuss the physical implications 
of our calculations, in particular in relation to different 
hadronic models of the nucleon. Results and conclusions 
are summarized in section 



II. CALCULATION OF THE CORRELATION 
FUNCTIONS FROM THE INSTANTON LIQUID 
MODELS. 

Our theoretical investigation is based on the analysis 
of two and three-point Euclidean correlation functions, in 
coordinate space. In particular we consider the electro- 
magnetic three-point function: 

Usc4{x,y) =< 0\tr[r]sc{x/2)J2"'{y)fisc{-x/2)-fi] |0 >, 

(1) 

where J|™(?;) is the fourth component of the electro- 
magnetic current, 

Jr(y) ^"(y)7M^(2/) -ld{y)j,diy). (2) 

risc{x) denotes a particular combination of loffe currents, 
that we refer to as the scalar current: 



risc{x) := f-abc[u°'{x) C75 d^{x)\ u^ix), 



(3) 



where a, 6, c are color indices. In principle, one can 
choose any linear combination of loffe currents, which 
excite states with the same addictive quantum numbers 
of the proton (as well as states with opposite parity) . Our 
choice is theoretically motivated by the fact that Tj'^ci^) 
contains explicitly a di-quark term, which was shown to 
receive maximal coupling to instanton's zero modes Q . 
The nucleon coupling to such current is defined as 



< 0\r]scix)\p, s >= Kscu{p, s)e" 



(4) 
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where s) is a Dirac spinor and s is the spin index. 

Let us speciahze now a particular choice]^ for the three- 
point function (|l|), in which: 



y ■ X — 0, and 
l:= \v\ = \x\. 



(5) 



We shall consider also the following two-point Eu- 
clidean proton correlation function 

n,,(x) 0|Tr[7/,,(x/2)77,,(-x/2)74] |0 >, (6) 

in which the propagation can been chosen along the Eu- 
clidean time axis. Due to the presence of 74, this cor- 
relator receives contribution from two quark zero-modes 
propagators in the field of an instanton ^. Notice that, 
in a theory with two massless flavors, this is the maximal 
number allowed by Pauli principle. 

Since the same factor A^^, which can not be deter- 
mined from experiments, appears both in the expression 
for the proton two-point function (|^), and in three-point 
hmction (flT), it is natural to consider the ratio: 



(7) 



in which such factor is divided out. In the next section 
we show that r(^) for I > 0.7 fm depends only on the 
the proton's mass and form-factors. At this scale, r(^) 
is dominated by non-perturbative effects and therefore 
represents a suitable tool to probe model descriptions of 
the non-perturbative strong dynamics. 

In instanton models, one assumes that the Euclidean 
QCD partition function is saturated by the configurations 
of an interacting ensemble of pseudo-particles (instantons 
and anti-instantons) which, in ordinary Minkowski space, 
are associated to the tunneling between degenerate clas- 
sical vacua. 

Each configuration of the ensemble is determined by 
a set of collective coordinates, specifying the the posi- 
tion, the size and the orientation in color space of each 
instanton. 

The parameters of the model are the average instanton 
size, p ~ 1/3 fm and the average instanton density, n ~ 
l/m~*. These values have been determined from several 
phenomenological studies , ||^ as well as from lattice 
simulations pl[ |. 

In the Random Instanton Liquid Model (RILM), one 
assumes that all instantons have the same size and that 
the distribution of all remaining collective coordinates is 



completely random. In the Interacting Instanton Liquid 
Model (IILM) the distribution is calculated from numer- 
ical simulation of the partition function, describing in- 
teraction between instantons. The most important and 
complicated part of this interaction is fcrmion exchanges 
between them, described by the fermionic determinant in 
the partition function. 

Both the RILM and IILM have been shown to have 
nucleon as a bound state, and give reasonable values 
for the its mass j^]. Furthermore, these models provide 
very different nucleon (octet) and delta (decuplet) corre- 
lators, explaining the splitting between the nucleon and 
the delta in terms of zero mode forces ||l2|, ||l3|. An- 
other important result of these works is the prediction 
of the nucleon coupling constants to loffe current: it is 
proportional to the probability to find all three quarks 
at the same point, and is needed for such applications as 
evaluation of proton decays. 

As the first step toward understanding of the function 
T{1) we have evaluated it analytically in the Single In- 
stanton Approximation (SIA), in which only the closest 
instanton is taken into account explicitly |2^], |^^. In 
such an approach, the correlation functions are evaluated 
from the the quark propagator S'/(^)(x, y) in the instan- 
ton (anti-instanton) background. This is known exactly 
to be the sum of two terms: a zero-mode part, Szmix, y) 
and a non-zero mode part, Snzm{x,y). In most appli- 
cations, Snzm{x,y) can be approximated with the free 
propagator, SQ{x,y), while the quark zero- mode propa- 
gator reads: 



Szm{x^ y) — 



i^o{x)ipl{y) 



(8) 



where iPq{x) represents an eigenmode of the Dirac oper- 
ator ^ with vanishing eigenvalue (quark zero-mode) , 



(9) 



ipo au{x; z) = 



7r((a;-z)2-f p2)3/2 

1-75 



(10) 



Q/3 



(here z denotes the instanton position, a, /3 = 1, 4 are 
spinor indices and Uab represents a general group ele- 
ment). 

In we showed that, if the correlation function re- 
ceives contribution from at least two quark zero-mode 
propagators, one can incorporate the effects of all other 



^Notice that, in such configuration, the vertices x/2, —x/2 
and y do not lie on the same line. The convenience of such 
choice will be discussed in the next section. 

^This can be verified by simply counting the quark chirality 
flips. 



■^The explicit expression for Tpo{x) is, of course, gauge de- 
pendent. The ILM is formulated in the singular gauge, where 
the topological charge is localized around each instanton. 



3 



instantons in an effective mass term that replacing the 
current mass in (^). More specifically, correlation func- 
tions receiving contribution from two zero-mode propaga- 
tors will be proportional to inverse powers of the current 
quark mass, m^, which will be replaced by an effective 
parameter . Such parameter depends on the ensemble 



and has been determined to be 



(65 MeVy for the 



RILM and ml ~ (105 MeVf for the IILM. The resulting 
expressions for Ilsc{l) and n4sc(',0 are Q 



and 



n.c(0 



15 12 np^ 



[g^ + {z^-lY + p2]3[^ + zl + p^]3' 



(11) 



96 



25 7r8;12 ^6 15^8^2 



d^z 



[{1/2 + zif + z^+ p^fl^ [{1/2 ~ z^f + z^ + p 



2 I ^2 I ^2l3/2 



corresponding ILM results up to quite large distances 
0.7-^ ~ \fm, depending on the correlation function) 
so there may be a window for such an approach in the 
present calculation. In figs. H and |l^ we have plotted the 
ratios Iisc{l) /^sc free{l) and n4sc(^0/n4 sc free 

(/, /) ob- 
tained in the SIA and from numerical simulations in the 
RILM and IILM |^. It is quite surprising to notice that 
the SIA results for three-point functions reproduce the 
corresponding ILM result all the way up to 1.2/to. For 
larger distances statistical errors on ILM results make 
the comparison somewhat meaningless. 

Now we make the second step, and numerically eval- 
uate the function T{1) from a random and an interact- 
ing instanton ensembles. In both cases we used propa- 
gators which include complete re-diagonalization in the 
zero-mode subspace. We have used 256 instantons and 
anti-instantons in a (4/m)'' periodic box. The results are 
presented in fig. |^: one can see that, although the IILM 
and the RILM give different two-point and three-point 
correlation functions (see figs, ^and |l^), their ratio r(/) 
is essentially the same. 



[{1/2 + Z4)2 + + p2]3/2 [(;/2 - Z^f + Z^ + p'^f'^ 

7/5 

[{I - Z^Y +Z^+ p2]3/2 _ ^2)2 + ^2]l/2 

f 2P + 5lz4 + + 2zl-Alz2 

' V [(^2 - + ^ + P^f^ [{1/2 + ZiY + ^ 

2P-5lz4 + ^ + 2zl-4:lz2 \] 
[{1/2 + z^y + + p2]3/2 [(;/2 _ ^4)2 + ^]l/2 ] J' 



(12) 

where := z2 + z| + z| and := z| + z| 4- zf. 

Note, that contribution from all other instantons en- 
ters only via an effective mass parameter. Although the 
values of these effective masses in RILM and IILM are 
quite different. Fortunately, one finds that two and three- 
point functions have the same power of this parameter, 
and thus it actually drops out in our ratio T{1) as soon as 
the perturbative (free) contribution becomes negligible. 

The SIA is supposed to break down when quarks 
are propagating for distances larger than the typical 
separation between two neighbor instantons (~ Ifm). 
Nevertheless, in |2^ we showed that the SIA calcula- 
tions for hadronic two-point functions agree with the 
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FIG. 3. Large distance behavior of T{1). The lines rep- 
resent the contribution of the nucleon pole, for the different 
scenarios given by (^o|) and (^l|). In particular, the solid line 
corresponds to the scenario 1, where the ratio hGe/Gm sta- 
bilizes at around 0.6 at high Q^, while the large-dashed line 
was obtained using the dipole formula for both Gm and Ge-- 
The circles (squares) represent RILM (IILM) predictions. 



*Our theoretical calculations for n4sc{l,l) have been ob- 
tained retaining only connected diagrams. This is motivated 
by the fact that all one-instanton contributions to the discon- 
nected diagrams for the electro-magnetic three-point function 
are mqzero. 



^In producing these plots we have used the effective mass 
parameters defined and evaluated in p3|. 
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III. RELATING THE EUCLIDEAN 
CORRELATION FUNCTIONS WITH THE 
EXPERIMENTAL DATA 

In a field theoretic Euclidean framework (such as lat- 
tice QCD or instanton-based models) a very important 
step is to encode available information about masses, 
form-factors, wave functions and other properties of par- 
ticular hadrons into appropriate Green functions. 

The spectral representation of correlators can be 
traditionally written assuming a "pole-plus-continuum" 
ansatz for the spectral density: 



1 1"°° 



(13) 

Here the first term represent the nucleon contribution: 
A/jv is the nucleon mass and D{m, x) is the Euclidean 
scalar propagator of a particle of mass m. 

The second term represents the contribution of the ex- 
cited states with nucleon quantum numbers (including 
those with opposite parity) . This particular form of "con- 
tinuum contribution" represents parton-hadron duality, 
it corresponds to the discontinuity of the perturbative 
(free) quark loop diagram contributing to the same cor- 
relation function. The new parameter Sq identifies the 
threshold for the onset of the duality regime: at this 
energy transition from the non-perturbative to the per- 
turbative regime takes place. The QCD sum rules practi- 
tioners obtain the value for this parameter together with 
others (masses and coupling constants), from comparing 
theoretical and experimental dependence of the correla- 
tor. However, the results still depend on this specific 
assumption about the contribution of the excited states. 

We have chosen not to deal with this uncertainty, cal- 
culating all correlation functions at sufficiently large dis- 
tances X <^ where the continuum contribution is 
negligible compared to nucleon contribution, (see fig. |^). 

The contribution to (|l|) due to the nucleon pole can be 
evaluated in terms of the Fourier transforms of the Dirac 
and Pauli form-factors (see appendix |^ for the detailed 
derivation) : 



Ilsc4:yx,y) ^ I d zTr 



S 



/x 



(14) 

where Sm,^{x,x') denotes the nucleon propagator from 



to X and 



F!^{z) ■.= ^^"F2{z)^—F^{z) 



(15) 



Eq. (|T^ has a simple physical interpretation (see fig. 
^): at large distances, the three-point function (|l|) rep- 
resents the (Euclidean) amplitude for a finite-sized nu- 
cleon to be created at —x/2, to absorb a photon at y and 



be annihilated at x/2. Notice that the electro-magnetic 
vertex depends on the Fourier transforms of the Pauli 
and Dirac form-factors and therefore contains informa- 
tion about the shape of the proton. 




FIG. 4. In the large distance limit, the three-point func- 
tion HscA is dominated by the nucleon pole. 

In order to evaluate Wsci{x,y)^ we need an expres- 
sion for Fi{z) and F2{z). This can be obtained by 
Fourier transforming the analytic continuation to Eu- 
clidean space of the fits of the experimental form-factors. 
Typically, elastic scattering experiments allow to obtain 
the Sachs electric and magnetic form-factors, which are 
related to the Pauli and Dirac form-factors by: 



G£(g2) =Fi(Q2)-T^^2(Q') 

Gm(Q') -i^i(0') + i^2(Q'), 



(16) 
(17) 



where r := Q'^/ANP and ^ _g2^ 

At low momenta, < i GeV"^, Ge{Q'^) and Ga/ (Q^) 
can be extracted from usual Rosenbluth separation 
method and error-bars are of the order of a few percent. 
In this regime, the data for both form-factors are well 
fitted by the so-called dipolc formula, 

GM(Q^)/f^^GE{Q^)^GD{Q^), Q^<lGeV^ (18) 



where 



Gz,(g2 



(1 + 1^)^^ 



(19) 



and /i is the proton magnetic moment. 

As the momentum increases, the cross section becomes 
dominated by Gm making the extraction of Gb more and 
more difficult. As a result, error-bars on G^ become very 
large already at ~ 2 GeV^ while error-bars on Gm 
remain smah up to ^ SQGeV^ |l|]. 

Using this experimental information we are now able to 
complete our phenomenological estimate of the relevant 
three-point function. The coordinate representation of 
Gm has been obtained by Fourier transforming the dipole 
fit (|l^), which is very accurate at the momentum scale 
we are interested in. The fit for Ge can be obtained as 
follows. One can assume the following ansatz: 



GE{Q')^GD{Q')aQ^), 



(20) 
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where the function ^(Q^), which accounts from the devi- 
ation from dipole shape, can be obtained from the avail- 
able data for /X /Gm ■ These can be well reproduced by 
the following set of exponential fits which, however, have 
different trends in the kinematical region > 3.5 GeV^ 
( see fig. |l| ): 

S.diQ'^) — 1- (dipole parametrization) 



= 0.39e 



+0.61. 

(q2)1.3 



^2(Q^) = 0.70 e — 53(5- + 0.30 



(21) 



The corresponding functions Fi{z) and ^2(2) are plotted 
in fig. 1^ and fig. ^. 




FIG. 5. Fits of the Pauli form-factor {F\) in coordinate 
representation obtained using the parametrization (Ell) 




distances, where the proton pole is dominant. We there- 
fore need to estimate how large I needs to be for this 
assumption to hold. From the analysis of the two-point 
function, we have seen that the continuum contribution 
is almost completely suppressed for I ^ 1 fm. In the 
electro-magnetic three-point function, however, one of 
the quarks is struck in the point y, acquiring and un- 
determined amount of momentum. As a result, at such 
vertex the nucleon pole gets mixed with all other res- 
onances and one needs to wait for another ~ 1 /m, in 
order for the pole to decouple again. From this rough es- 
timate, one expects that (||) can be identified with ( p^ if 
the three quarks have propagated for about 2/m. In the 
ILM and on the lattice, the vacuum is usually simulated 
in a four-dimensional periodic box, whose maximal size 
is determined by the machine's precision and by com- 
putational time constraints. It is therefore important to 
make a convenient choice of x and y in (|l]). In partic- 
ular, we need to consider three-point functions that are 
large enough for the ground-state pole to decouple but 
can be fit into a relatively small box. At this purpose, it 
is generally a good idea to let quarks propagate along the 
diagonals. With the choice (p) one has that the quarks 
travel for the needed ^ "ifm already for I ^ 0.7 fm. 

At this point, we may ask to what momentum scale is 
our analysis in configuration space sensitive to. From un- 
certainty principle arguments, one would naively expect 
that correlation functions of size I should be essentially 
sensitive to the physics at the scale Q ~ l/l. Along 
this line, one is then lead to argue that our large sized 
{I > 0.7 fm) correlation functions can only be used to 
study to the small momentum part of the form-factors. 
This is certainly true asymptotically: for example, in 
the (numerically infeasible) limit of extremely large /, 
the three-point function should depend only on the pro- 
ton electric charge and magnetic moment and, therefore, 
all scenarios considered in (|2l|) must give the same an- 
swer. However, at intermediate distances the situation is 
more complicated. In fact, eq. (14) shows that nsc4(0j 
is related through a non-trivial i nteg ral to the Fourier 
transform of the form- factors (eq. |l4| ). Such integration 
mixes the different modes in such a way that three-point 
functions with I > 0.7 fm are still quite sensitive to the 
shape of form- factors at high Q^. This is already seen in 
fig. 0, where Il4sc{l,l), obtained form fits of the form- 
factors that differ only in the high momentum range ( 
> 3.5 GeV^) are compared. Clearly, such green func- 
tions at, say, / ~ 1.5 fm are quite different and larger 
distances are needed for the three curves to converge. 



FIG. 6. Fits of the Dirac form-factor {F2) in coordinate 
representation obtained using the parametrization ( pl| ) 



The results of our phenomenological calculation for 
T{1), using the parametrization ( pi| ) for £,{Q^) are re- 
ported in fig. m It is important to recall that our phe- 
nomenological calculation of T{1) is valid only at large 
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between various "scenarios". In fact, the scenario (in 
which GE^J■/GM ~ 1, identically) and the scenario 3 
(in which Ge^J-IGm rapidly approaches zero) seem to 
be ruled-out, assuming the theory prediction is correct. 
Best agreement seem to be reached for the scenario 1, 
the case in which Ge IJ^/Gm falls with till about 0.6, 
as found by Jlab experiment, and then more or less sta- 
bilizes at large momenta. Scenarios in which the ratio 
stabilizes at similar values are of course also possible. 



FIG. 7. Proton pole contribution to the Electro-magnetic 
three-point function 114 ac, obtained from three different fits 



for ^ Ge I 



21h which differ only for Q > 3.5 GeV^ 



It is interesting to further analyze the contribution of 
different momentum components to our correlator. At 
this purpose, in fig. ^ we have plotted the three-point 
function corresponding to scenario 1, and we compare it 
with two more curves, that one obtains from the same 
parametrization, but keeping only the contribution from 
Q < IGeV and Q > IGeV to the form-factors Fi{z) 
and ^2(2) (see eq. A7). From this figure we clearly see 
that some non trivial cancellations occur and the high 
momentum components of the form factors play an im- 
portant role in the whole range of I considered. 
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FIG. 9. Proton two-point function, Ilsc{l) evaluated in 
the ILM and in the SIA. The circles (squares) correspond to 
RILM (IILM) results. Dashed (dot-dashed) lines correspond 
to SIA results using the mass parameters evaluated in [ p3[ : 
mi = {65 MeV f ( ml = {107 MeV f ). 
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FIG. 8. Contribution of different momenta to the proton 
electro- magnetic three-point function 114 sc ■ The solid line 
represents the correlators obtained in scenario 1, integrating 
over all momenta, the dotted (dashed) line is obtained keep- 
ing only the Q > 1 GeV {Q < 1 GeV^) components of the 
form-factors. 

Now we are ready to compare theoretical calculations 
for (^ to the four curves corresponding to different "sce- 
narios" for the electric form- factor behavior at large . 
First of all, one finds that the overall magnitude and 
qualitative behavior with I are similar. This is very non- 
trivial, since there are compensations between electric 
versus magnetic, and small versus large . 

Furthermore, our results seem to be able to distinguish 



0.2 0.4 0.6 0.8 I 1.2 

I [tm] 

FIG. 10. Proton three-point function, Il4,sc{l) evaluated 
in the ILM and in the SIA. The circles (squares) correspond to 
RILM (IILM) results. Dashed (dot-dashed) lines correspond 
to SIA results using the mass parameters evaluated in [ p3[ : 
ml = {65 MeV f { ml = {107 MeV f ). 



IV. DISCUSSION AND COMPARISON WITH 
HADRONIC MODELS 

After having calculated of the relevant Green func- 
tions, analytically in SIA and numerically for the ensem- 
bles, and compared our prediction with the data, let us 
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go into somewhat more general discussion of the physics 
involved. 

First let us ask, why RILM and and IILM give so close 
results for the ratio T{1), but not the correlators them- 
selves? The former model has no correlation between 
instantons while the latter has strong ones^ and their 
predictions are usually not that close. The reason is that 
(see also discussion in section ||), at the scale we are 
working at, the correlation functions are dominated by 
one single instanton. 

This implies two important consequences. First, there 
actually exists a "working window" where the (unknown) 
continuum contribution is highly suppressed and yet SIA 
is still accurate. Second, since the only relevant dimen- 
sionful parameter in SIA is the instanton size p (the in- 
stanton density also drops out in the ratio), it follows 
that the scale of the nucleon form-factor is completely 
determined by itj^. 

The content of the typical nucleon configuration is pic- 
tured phenomenologically in different ways: e.g. the 
usual quark model picture is that of three massive "con- 
stituent" quarks, bound by confining forces. Another 
possibility discussed in literature is as a loose bound state 
of a quark with a scalar system of two quarks, more 
tightly bound (di-quark). In this section, we consider 
the problem of qualitatively determining the size of such 
di-quark, in the instanton-based models. 

From the Pauli principle on the level of instanton zero 
modes, only two (u and d) of the three quarks can be in 
the state described by the zero-mode wave-function ( |l0| ) 
and therefore be localized around the closest instanton. 
As a result, the typical size of the di-quark should be 
determined by the size of the instanton, and be generally 
comparable to it. 

There is an interesting theoretical argument suggested 
in |2j], which relates pions and scalar di-quarks in the 
QCD with two colors. Indeed, in this theory di-quarks 
are not much heavier than mesons, but are degener- 
ate with them, due to the so called Pauli-Gursey sym- 
metry. In particular, the lowest di-quarks should be 
bound as strongly as the lowest meson, the massless pi- 
ons. And indeed, the general pattern of symmetry break- 
ing SU{2Nf) Sp{2Nf) and the number of Goldstone 
modes is 

Ngoldstones = "^Nj — Nf — 1. 

For the most interesting case Nf = 2 there are five mass- 
less modes: three pions, plus scalar di-quark S and its 



anti-particle S. So, in this theory the pions and scalar 
di-quarks are related by symmetry, and should be truly 
identical. For three colors it is of course not so, but by 
continuity in the number of colors, one may expect some 
traces of that to remain valid. 

To get quantitative estimates, we proceed as follows. 
On the one hand, we evaluate from the ILM the di-quark 
two-point function, 

n,,(0=<0| J,'',(x)JL,(0)|0>, (22) 
where J^ci^) is the scalar current: 

J^M^e'^''^u\x){C-f5)d%x), (23) 

and the di-quark three-point function, 

U,,,{l,l) Omx/2) Jriy) JLA^x/2)\0 >, (24) 

and obtain the ratio Tdqil) := Il4dq{l'J)/^dq{l)- 

On the other hand, we compute the spectral decom- 
position for the same quantities, as if the di-quark be a 
physical (scalar) particle. This assumption was discussed 
in detail in , where the mass of the di-quark was also 
determined from the RILM to be Mdq = 420 ± 30 MeV. 
As in the case of the nucleon, the large distance limit 
of the spectral decomposition of the three-point function 
requires the knowledge of the form-factor of lowest ly- 
ing state. Exploiting the similarities between the phys- 
ical properties of the pion and those of the di-quark in 
the ILM, we take a mono-pole shape ansatz for its form- 
factor, 

FdqiQ') = (25) 

The free parameter a can be determined by comparing 
the function Tdq evaluated (at large distances) from ILM 
and form the spectral decomposition (see fig. |ll|) we 
found a - (ITOOMeV)^. 

Clearly, the agreement of ILM prediction with the pole 
contribution to the spectral representation (obtained us- 
ing (psf)) is definitely worse for the scalar di-quark than 
for the pion |p^ , indicating that the simple mono-pole 
ansatz ( p5| ) does not quite work in the case of the di- 
quark. Of course this should not surprise: after all, the 
di-quark is not an asymptotic state and we expect its 
charge distribution to be somewhat distorted by the pres- 
ence of the third quark. 



To remind the reader about importance of these correlation: 
in IILM topological charge is screened and topological sus- 
ceptibility is small, 0{mq), as in the QCD vacuum with light 
quarks. In RILM it is large, like in gluodynamics (quenched 
QCD). 

^Note that this feature is similar to what has been found for 
the pion form-factor in jl6|. 
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FIG. 11. The ratio Tsc4{l,l) evaluated in the ILM, SIA 
(solid line) and from the ansatz ( |25[ ) with a = 1760MeV^ 
(dashed line). 

Nevertheless, from this fit we can still obtain a quali- 
tative estimate of its size from the usual relation: 



< >:= 



6 



0.3/m. 



(26) 



This value can, for instance, be compared with the phe- 
nomenological quark-diquark model | ^ , where the value 
< r^^ >~ (0.45 /to)^ was found. 

Although both values can only be considered as qual- 
itative estimates, they seem to indicate that the size of 
the di-quark is of the same order of the typical size of an 
instanton. 



We have also performed a similar analysis to obtain 
a qualitative estimate of the size of the di-quark part 
of the nucleon. We found < (0.3 /m)^, which is 

consistent with the picture in which quark-zero modes are 
localized around the closest instanton and in reasonable 
agreement with the value < >~ (0.45 fm)^ found in 

Our analysis revealed that, although the RILM and 
the IILM give in general quite different correlators, their 
prediction for the ratio T{1) is exactly the same. This 
implies that at the moment of scattering two out of three 
quarks interact essentially with one instanton, via its zero 
mode. This particular configuration is, in a way, at a core 
of quark-diquark model. 

As we found that it is quite possible to perform an anal- 
ysis of the pion and proton form-factor in the simplest 
SIA, without needing to model the continuum contribu- 
tion, one can hope that such analytical results can be 
Fourier transformed to momentum space and allow for 
direct comparison with experimental data. We intend to 
exploit this possibility in our further works. 
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V. CONCLUSIONS AND OUTLOOK 



In this paper we have studied the instanton contribu- 
tion to the proton's electro-magnetic form- factors. We 
have calculated the ratio TscQ) of the proton electro- 
magnetic three-point function (|l|) and the two-point func- 
tion (^) in the RILM, in the IILM and, analytically, in 
the SIA. It is important to stress that there is no free 
parameters and therefore the results of fig. |^ represent 
an absolute prediction for this Green function. Further- 
more, the results happen to be independent on the de- 
tails of the model, and (like it was previously found for 
the pion form-factor) only the instanton mean radius 
p ^ 1/600 Mew is actually setting the scale of it. 

We have compared our theoretical predictions with the 
phenomenological curve obtained from the nucleon's pole 
contribution to spectral representation of the correlation 
functions. In particular, in order to evaluate the proton 
contribution to the electro-magnetic three-point function 
we have used several fits of the Ge^J^/Gm data. These 
fits describe scenarios in which such ratio reaches differ- 
ent asymptotic values, in the kinematical region under 
current investigation at JLAB. We found that our pre- 
diction, is consistent with the deviation from the dipole 
shape of and is very much in agreement with a sce- 
nario in which the Ge p/Gm ratio stabilizes somewhere 
around 0.6, at high Q"^ . 
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APPENDIX A: DERIVATION OF THE 
NUCLEON POLE CONTRIBUTION TO THE 
ELECTRO-MAGNETIC THREE-POINT 
FUNCTION 

We want to evaluate the proton's pole contribution to 
the three-point function: 

U,,i{x,y) 0\Tr[rj,,ix)jr{y)fisMli] |0 >, (Al) 
where J|™ is the electro-magnetic current: 

jriy) ■■= I uiy) 74 uiy) - i J(y) jM, (A2) 

. and we haven chosen X4 > > 0. 

We begin by inserting two complete sets of states: 



<rV4(y)|r><r|,7,c(-^) |o>74 



(A3) 



The term of this double sum that depends only on the 
lowest lying state (nuclcon) represents the pole, while the 
sum of all other terms forms the continuum. We therefore 
obtain: 



p 



< 



2wp, (27r)3 J 2w;,(27r)3 

■7^H<0hsc (I) \p\s';Mn >■ 
■ < p', s'; Mn I My) I P, s; Mn > ■ 
p,s;A/Ar|77,c(-|) |0 > 74] + --, 



(A4) 



where Mn denotes the proton's rest mass and the ellipses 
denote the continuum contribution. Now, using equation 
and recalling the definition of Pauli and Dirac form- 
factors, 

< p', s'; Mwl J^(0)|p, s;Mn> = 
= u{p\s') [7^Fi(Q2) + l^a^,q'^F2{Q^)]u{p,s), (A5) 



we get: 



'^sc.i{x,y) 



d'p 



2ijjp{2iTf J 2tj^, (27r) 



exp 



'^I^{p + P')~^y{p~p')\ Tr[{P' + MN)- 



(A6) 



Now we can re-write the form-factors in terms of their 
Fourier transform, 



d^z 



(A7) 
(A8) 



and introduce the expression for the nucleon propagator 
in coordinate space {x'q>xq ): 



obtaining: 



d^p 



(27r)3 2tjp 



-^P-(^'-^){f)+MN), (A9) 



n,c4(2;,y) = Af, / d^zTr 



S 



Sm^ (y + z,-^^ 74 



(AlO) 



where: 



F^iz) ■.= td''F2{z). 



(All) 



One can check that the analytic continuation of ( A10| ) 
to Euclidean space ^ is real. 
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